Error analysis for mesospheric temperature profiling by absorptive occultation sensors

International audienceAn error analysis for mesospheric profiles retrieved from absorptive occultation data has been performed, starting with realistic error assumptions as would apply to intensity data collected by available high-precision UV photodiode sensors. Propagation of statistical errors was investigated through the complete retrieval chain from measured intensity profiles to atmospheric density, pressure, and temperature profiles. We assumed unbiased errors as the occultation method is essentially self-calibrating and straight-line propagation of occulted signals as we focus on heights of 50?100 km, where refractive bending of the sensed radiation is negligible. Throughout the analysis the errors were characterized at each retrieval step by their mean profile, their covariance matrix and their probability density function (pdf). This furnishes, compared to a variance-only estimation, a much improved insight into the error propagation mechanism. We applied the procedure to a baseline analysis of the performance of a recently proposed solar UV occultation sensor (SMAS ? Sun Monitor and Atmospheric Sounder) and provide, using a reasonable exponential atmospheric model as background, results on error standard deviations and error correlation functions of density, pressure, and temperature profiles. Two different sensor photodiode assumptions are discussed, respectively, diamond diodes (DD) with 0.03% and silicon diodes (SD) with 0.1% (unattenuated intensity) measurement noise at 10 Hz sampling rate. A factor-of-2 margin was applied to these noise values in order to roughly account for unmodeled cross section uncertainties. Within the entire height domain (50?100 km) we find temperature to be retrieved to better than 0.3 K (DD) / 1 K (SD) accuracy, respectively, at 2 km height resolution. The results indicate that absorptive occultations acquired by a SMAS-type sensor could provide mesospheric profiles of fundamental variables such as temperature with unprecedented accuracy and vertical resolution. A major part of the error analysis also applies to refractive (e.g., Global Navigation Satellite System based) occultations as well as to any temperature profile retrieval based on air density or major species density measurements (e.g., from Rayleigh lidar or falling sphere techniques)

Greenhouse gas profiling by infrared-laser and microwave occultation: retrieval algorithm and demonstration results from end-to-end simulations

Measuring greenhouse gas (GHG) profiles with global coverage and high accuracy and vertical resolution in the upper troposphere and lower stratosphere (UTLS) is key for improved monitoring of GHG concentrations in the free atmosphere. In this respect a new satellite mission concept adding an infrared-laser part to the already well studied microwave occultation technique exploits the joint propagation of infrared-laser and microwave signals between Low Earth Orbit (LEO) satellites. This synergetic combination, referred to as LEO-LEO microwave and infrared-laser occultation (LMIO) method, enables to retrieve thermodynamic profiles (pressure, temperature, humidity) and accurate altitude levels from the microwave signals and GHG profiles from the simultaneously measured infrared-laser signals. However, due to the novelty of the LMIO method, a retrieval algorithm for GHG profiling is not yet available. Here we introduce such an algorithm for retrieving GHGs from LEO-LEO infrared-laser occultation (LIO) data, applied as a second step after retrieving thermodynamic profiles from LEO-LEO microwave occultation (LMO) data. We thoroughly describe the LIO retrieval algorithm and unveil the synergy with the LMO-retrieved pressure, temperature, and altitude information. We furthermore demonstrate the effective independence of the GHG retrieval results from background (a priori) information in discussing demonstration results from LMIO end-to-end simulations for a representative set of GHG profiles, including carbon dioxide (CO<sub>2</sub>), water vapor (H<sub>2</sub>O), methane (CH<sub>4</sub>), and ozone (O<sub>3</sub>). The GHGs except for ozone are well retrieved throughout the UTLS, while ozone is well retrieved from about 10 km to 15 km upwards, since the ozone layer resides in the lower stratosphere. The GHG retrieval errors are generally smaller than 1% to 3% r.m.s., at a vertical resolution of about 1 km. The retrieved profiles also appear unbiased, which points to the climate benchmarking capability of the LMIO method. This performance, found here for clear-air atmospheric conditions, is unprecedented for vertical profiling of GHGs in the free atmosphere and encouraging for future LMIO implementation. Subsequent work will examine GHG retrievals in cloudy air, addressing retrieval performance when scanning through intermittent upper tropospheric cloudiness

Errors in GNSS radio occultation data: relevance of the measurement geometry and obliquity of profiles

Atmospheric profiles retrieved from GNSS (Global Navigation Satellite System) radio occultation (RO) measurements are increasingly used to validate other measurement data. For this purpose it is important to be aware of the characteristics of RO measurements. RO data are frequently compared with vertical reference profiles, but the RO method does not provide vertical scans through the atmosphere. The average elevation angle of the tangent point trajectory (which would be 90° for a vertical scan) is about 40° at altitudes above 70 km, decreasing to about 25° at 20 km and to less than 5° below 3 km. In an atmosphere with high horizontal variability we can thus expect noticeable representativeness errors if the retrieved profiles are compared with vertical reference profiles. We have performed an end-to-end simulation study using high-resolution analysis fields (T799L91) from the European Centre for Medium-Range Weather Forecasts (ECMWF) to simulate a representative ensemble of RO profiles via high-precision 3-D ray tracing. Thereby we focused on the dependence of systematic and random errors on the measurement geometry, specifically on the incidence angle of the RO measurement rays with respect to the orbit plane of the receiving satellite, also termed azimuth angle, which determines the obliquity of RO profiles. We analyzed by how much errors are reduced if the reference profile is not taken vertical at the mean tangent point but along the retrieved tangent point trajectory (TPT) of the RO profile. The exact TPT can only be determined by performing ray tracing, but our results confirm that the retrieved TPT – calculated from observed impact parameters – is a very good approximation to the "true" one. Systematic and random errors in RO data increase with increasing azimuth angle, less if the TPT is properly taken in to account, since the increasing obliquity of the RO profiles leads to an increasing sensitivity to departures from horizontal symmetry. Up to an azimuth angle of 30°, however, this effect is small, even if the RO profiles are assumed to be vertical. For applications requiring highest accuracy and precision it is advisable to exclude RO profiles with ray incidence angles beyond an azimuth of 50°. Errors in retrieved atmospheric profiles decrease significantly, by up to a factor of 2, if the RO data are exploited along the retrieved TPT. The tangent point trajectory of RO profiles should therefore be exploited whenever this is possible

Retrieval of temperature profiles from CHAMP for climate monitoring: intercomparison with Envisat MIPAS and GOMOS and different atmospheric analyses

International audienceThis study describes and evaluates a Global Navigation Satellite System (GNSS) radio occultation (RO) retrieval scheme particularly aimed at delivering bias-free atmospheric parameters for climate monitoring and research. The focus of the retrieval is on the sensible use of a priori information for careful high-altitude initialisation in order to maximise the usable altitude range. The RO retrieval scheme has been meanwhile applied to more than five years of data (September 2001 to present) from the German CHAllenging Minisatellite Payload for geoscientific research (CHAMP) satellite. In this study it was validated against various correlative datasets including the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and the Global Ozone Monitoring for Occultation of Stars (GOMOS) sensors on Envisat, five different atmospheric analyses, and the operational CHAMP retrieval product from GeoForschungsZentrum (GFZ) Potsdam. In the global mean within 10 to 30 km altitude we find that the present validation observationally constrains the potential RO temperature bias to be <0.2 K. Latitudinally resolved analyses show biases to be observationally constrained to <0.2?0.5 K up to 35 km in most cases, and up to 30 km in any case, even if severely biased (about 10 K or more) a priori information is used in the high altitude initialisation of the retrieval. No evidence is found for the 10?35 km altitude range of residual RO bias sources other than those potentially propagated downward from initialisation, indicating that the widely quoted RO promise of "unbiasedness and long-term stability due to intrinsic self-calibration" can indeed be realised given care in the data processing to strictly limit structural uncertainty. The results thus reinforce that adequate high-altitude initialisation is crucial for accurate stratospheric RO retrievals. The common method of initialising, at some altitude in the upper stratosphere, the hydrostatic integral with an upper boundary temperature or pressure value derived from meteorological analyses is prone to introduce biases from the upper boundary down to below 25 km. Also above 30 to 35 km, GNSS RO delivers a considerable amount of observed information up to around 40 km, which is particularly interesting for numerical weather prediction (NWP) systems, where direct assimilation of non-initialised observed RO bending angles (free of a priori) is thus the method of choice. The results underline the value of RO for climate applications

Small catchment runoff sensitivity to station density and spatial interpolation: Hydrological modeling of heavy rainfall using a dense rain Gauge network

Precipitation is the most important input to hydrological models, and its spatial variability can strongly influence modeled runoff. The highly dense station network WegenerNet (0.5 stations per km2) in southeastern Austria offers the opportunity to study the sensitivity of modeled runoff to precipitation input. We performed a large set of runoff simulations (WaSiM model) using 16 subnetworks with varying station densities and two interpolation schemes (inverse distance weighting, Thiessen polygons). Six representative heavy precipitation events were analyzed, placing a focus on small subcatchments (10–30 km2) and different event durations. We found that the modeling performance generally improved when the station density was increased up to a certain resolution: a mean nearest neighbor distance of around 6 km for long-duration events and about 2.5 km for short-duration events. However, this is not always true for small subcatchments. The sufficient station density is clearly dependent on the catchment area, event type, and station distribution. When the network is very dense (mean distance &lt; 1.7 km), any reasonable interpolation choice is suitable. Overall, the station density is much more important than the interpolation scheme. Our findings highlight the need to study extreme precipitation characteristics in combination with runoff modeling to decompose precipitation uncertainties more comprehensively

Refractivity and temperature climate records from multiple radio occultation satellites consistent within 0.05%

Data consistency is an important prerequisite to build radio occultation (RO) climatologies based on a combined record of data from different satellites. The presence of multiple RO receiving satellites in orbit over the same time period allows for testing this consistency. We used RO data from CHAMP (CHAllenging Minisatellite Payload for geoscientific research), six FORMOSAT-3/COSMIC satellites (Formosa Satellite Mission 3/Constellation Observing System for Meteorology, Ionosphere and Climate, F3C), and GRACE-A (Gravity Recovery and Climate Experiment). We show latitude-altitude-resolved results for an example month (October 2007) and the temporal evolution of differences in a climate record of global and monthly means from January 2007 to December 2009. Latitude- and altitude-resolved refractivity and dry temperature climatologies clearly show the influence of different sampling characteristics; monthly mean deviations from the multi-satellite mean over the altitude domain 10 km to 30 km typically reach 0.1% and 0.2 K, respectively. Nevertheless, the 3-yr average deviations (shorter for CHAMP) are less than 0.03% and 0.05 K, respectively. We find no indications for instrument degradation, temporal inhomogeneities in the RO records, or temporal trends in sampling patterns. Based on analysis fields from ECMWF (European Centre for Medium-Range Weather Forecasts), we can estimate &amp;ndash; and subtract &amp;ndash; the sampling error from each monthly climatology. After such subtraction, refractivity deviations are found reduced to &lt;0.05% in almost any month and dry temperature deviations to &lt;0.05 K (&lt;0.02% relative) for almost every satellite and month. 3-yr average deviations are even reduced to &lt;0.01% and &lt;0.01 K (CHAMP: &amp;minus;0.05 K), respectively, establishing an amazing consistency of RO climatologies from different satellites. If applying the same processing scheme for all data, refractivity and dry temperature records from individual satellites with similar bending angle noise can be safely combined up to 30 km altitude (refractivity also up to 35 km) to a consistent single climate record of substantial value for climate monitoring in the upper troposphere and lower stratosphere

Where to see climate change best in radio occultation variables – study using GCMs and ECMWF reanalyses

Radio occultation (RO) is a new technique to observe the upper troposphere and lower stratosphere (UTLS), a region that reacts particularly sensitive to climate change. Featuring characteristics such as long-term stability, SI traceability, all-weather capability, global coverage, and high accuracy and vertical resolution, RO data fulfill the requirements for climate monitoring in the UTLS. However, while a range of studies has shown the climate utility of RO it has not yet been explored sytematically where to see climate change best in RO variables. Therefore we perform here a systematic trend study for the RO variables refractivity, pressure, and temperature (bending angle, not depending on height but impact parameter, is left for separate study). The trends, given at geopotential height levels and for layer gradients, are explored to determine seasons, geographic regions, and height domains, which show a significant trend signal. Because continuous RO data are available since 2001 only, reanalyses (ERA-40 and ERA-Interim) and global circulation model simulations of the Intergovernmental Panel on Climate Change Assessment Report 4 (CCSM3, ECHAM5, HadCM3) are used as proxy data for RO. It is shown that RO data are sensitive at different height ranges and that thus several indicators of climate change can be retrieved. Refractivity emerges as indicator in the lower stratosphere (LS) and tropopause region at about 14 km to 24 km, pressure over the whole UTLS, and both in all large-scale regions except the polar caps. Temperature qualifies as indicator in the upper troposphere below about 16 km and in the LS above about 21 km. Overall, refractivity and pressure alone are adequate indicators for the UTLS, but temperature as commonly used variable facilitates easy interpretation of results. Layer gradients were found to be further sensitive indicators providing additional information. Besides large-scale global and hemispheric means the tropics and the mid-latitudes appear as regions suitable to track climate change with RO data. The results also point to the value of utilizing in addition to annual means specific seasons, such as northern hemispheric fall and summer, for early climate signal detection. Since RO data feature much better vertical resolution than the proxy data of this study, more detailed insights can be expected when a longer RO record will be available